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II. AMENDMENTS TO THE CLAIMS : 

Please cancel claims 2 and 3 without prejudice. Kindly amend claims 4-12, 14 and 15 
as follows. 

The following Listing of Claims replaces all prior listing, or versions 
Listing of Claims: 

1 . (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, directly using V-CAD data, comprising the steps of : 



boundaries orthogonal to each other, the external data including boundary data of an object 
that whieh contacts incompressible viscous fluid; 

a c e ll cl assifying s tep (B)-ef classifying the divided cells into an internal cell (13a) 
positioned inside or outside the object and a boundary cell (13b) including the boundary data; 

a cut point determining step (C)-ef determining cut points in ridges of the boundary 
cell on the basis of the boundary data; 

a bounda ry face de t e rminin g st e p (D) e£determining a polygon connecting the cut 
points to be cell internal data for the boundary face; and 

a analyzing step (E) ef-applying a cut cell finite volume method combined with a 
VOF method to a boundary of a flow field to analyze the flow field , wherein step (E) 
comprises the steps of 

i. applying a two-dimensional QUICK interpolation scheme to a convection 
term for space integral; 

ii. applying a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combining the convection term and the diffusion term, and applying an 




(A)-ef dividing external data into a plurality of cells (13) having 
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Adams-Bashforth method having precision of a degree of a second order to the 
combined convection term and diffusion term for time marching; and 

iv. applying a Euler implicit method having precision of a degree of a first 
order to a pressure gradient term for time marching, 
wherein for a two-dimensional boundary cell, a governing equation in the finite 
volume method is expressed by a governing equation (7), 



jj^F = - jjdiv(u ®u)dV- \\div{pl )dV + — \\div{grad{u))dV ( 7 ) 



wherein v *- f corresponds to the convention term, 



corresponds to the pressure gradient term, and 



1 

Re 



corresponds to the diffusion term . 



2. (Cancelled) 

3. (Cancelled) 



4. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 claim 3 , wherein the convection term, the 
pressure gradient term and the diffusion term in the governing equation of the finite volume 
method are expressed by the equations (8), (9) and (10) of Formula 2 , respectively, 
[Formula 2] 
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convection term: 

jfdiv(u®u)dV = <§(u®u) ndS = £ (Q®u) m n8S m 

y.. S,_ 5 "=1-5 

=[AMA^y-.y-^-.x:> I -w) 

+Ax(A, J u}» nj , V2 v, J - A^uX^v,^ )]7 ( 8 ) 

+[AX5,,v (x) « i ,-5,. 1 ,v w «,._,,) 

«- ^ v ,+1/2,/-]/ 2 i-\UJ-\(2 

pressure gradient term: 

\\div{pI)dV = <§{plyndS= £ A.' ^ 

- &y{B L jP M/2J - B„ jPt _ U2J - Pp (*,, - S,._ Iiy . )]/ ( 9 ) 

+ M4,A-, + i/2 -Aj-iPij-m-PpMj -4./-i)I7 
diffusion term: 

jjdiv(grad(u))dV = (|J grad{u) • rttfS = J] g™d(u) m • «^5 OT 

=[ak*,.,^(^ -(*,,. -5 M . y )^(i/)p 

(10) 

+Ax(4 jg rad(vX J+U2 - ^gradiyY^ -(4 . -4. y -,)^(v);)]7 

5. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 claim 3, wherein when a no-slip boundary 
condition is used for a solid boundary, a first integral is performed on the solid boundary with 
the convection term being zero, a value of a middle point P of a cut line segment isbemg used 
as an average for the pressure gradient term and the diffusion term, and a space integral is 
performed with areas fractions being applied to all of the terms. 



6. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 claim 3 , wherein the boundary cell having 
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the parameter smaller than a threshold value of VOF=0.01 is regarded as a complete solid, 

for the boundary cell having the parameter larger than the threshold value, a 

definition point for the parameter calculated in a cut cell is set at a center of the boundary 

cell, 

and a definition point for a parameter in a ridge is set at a center of a cell ridge, and a 
parameter at a middle point of a line segment-4 is calculated by a linear interpolation. 

7. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 claim 3 , wherein a drag force (in a flow 
direction) and a lift force (in a direction vertical to the flow^) acting on the object, are 
expressed by equations (12) and (13) of Formula 3, that respectively express drag force and 
lift force as follows. 



[For mu la 3] 
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drag force: 

= jj(^. )dxdy+ jjc^ydydc^a^ds+cj^cb (12) 

y ox y oy s s 

= JK* \f xW \f 2 iy)My+ J(^y \ SlM ^^{only Cartesian 

lift force: 

v v y s s 



8. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 1 claim 3 , wherein in fluid-structure 
interaction analysis accompanying a moving boundary, a fluid system and a structure system 
are separately analyzed for each predetermined time interval, and boundary conditions for the 
fluid system and the structure system are explicitly used. 



9. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 8, wherein in analysis on a forcibly vibrated 
circular cylinder, the circular cylinder is set as a_one-mass-point and a^one-degree-of- freedom 
system such that the circular cylinder is a solid structure elastically supported and vibrating in 
a direction vertical to the flow, 

and Y-direction displacement of a center of the circular cylinder is given by4fee 
equation (17), and a velocity boundary condition in the Y direction for a surface of the 
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circular cylinder is given by the-equation (18), 
y=Asin(2 7T f c t) ...(17) 
v w =A2 it f c cos(2 7i f c t) ...(18). 

10. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 9, wherein movement velocity of the 
vibrating circular cylinder obtained by the equation (18) is changed to be given for each 
calculation time step for the velocity boundary condition on the flow field. 

1 1 . (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 8, wherein in analysis on self-induced 
vibration due to aan vortex shedding from the circular cylinder, a vibration equation having a 
dimension is expressed by fee-equation (19) or equation (20), using a^one-mass-point and a 
one-degree-of-freedom dumper/spring model, 

m ^ +c ^ +k y={p u « DC L < 19 > 

^^¥J^H2,fSy~c L (20) 

12. (Currently Amended) A method for numerical analysis of a flow field of 
incompressible viscous fluid, according to claim 11, wherein the movement velocity of the 
vibrating circular cylinder calculated by tha-equation (20) is changed to be given for each 
calculation time step for the velocity boundary condition on the flow field. 

13. (Original) A method for numerical analysis of a flow field of incompressible 
viscous fluid, according to claim 11, wherein initial displacement and initial velocity of the 
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circular cylinder are set to be zero, the lift force is explicitly given by using a current value, 
and the vibration equation is integral by the Newmark's j3 method to obtain vibration 
displacement and vibration velocity of the circular cylinder. 

14. (Currently Amended) A device for numerical analysis of a flow field of 
incompressible viscous fluid, directly using V-CAD data, the device comprising: 

an input device~(2) for inputting external data including boundary data of an object-(4) 
that contacts incompressible viscous fluid; 

an external storage device-(3) for storing substantial data of shape data and physical 
property data integrated into each other, and a storage operational program for the substantial 
data; 

an internal storage device-(4) and a_central processing device-^ for executing the 
storage operational program; and 

an output device-{6) for outputting a result of fee-execution of the storage operational 
program; 

wherein the devic e for numerical analysis 

(A) divides the external data into a plurality of cells (13) having boundaries 
orthogonal to each others 

(B) classifies the divided cells into an internal cell (13a) positioned inside or 
outside the object and a boundary cell (13b) including the boundary data; 7 

(C) determines cut points in ridges of the boundary cell on the basis of the 
boundary data; ? 

(D) determines a polygon connecting the cut points to be cell internal data for 
the boundary face^ and 

(E) applies a cut cell finite volume method combined with a VOF method to a 
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boundary of a flow field to analyze the flow field , wherein when the device for 

numerical analysis applies the cut cell finite volume method combined with a VOF 

method to the boundary of the flow field to analyze the flow field, the device for 

numerical analysis operates to 

i. apply a two-dimensional QUICK interpolation scheme to a 
convection term for space integral; 

ii. apply a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combines the convection term and the diffusion term, and applies 
an Adams-Bashforth method having precision of a degree of a second order to 
the combined convection term and diffusion term for time marching; and 

iv. applies a Euler implicit method having precision of a degree of a 
first order to a pressure gradient term for time marching, 
wherein for a two-dimensional boundary cell, a governing equation in the finite 
volume method is expressed by a governing equation (7), 

jj^ictV = - jjdiv(H ® u W - jjdiv(pI)dV + — \\div(grad(u ))dV ( 7 ) 

— i 

-jfdivis mytv - lldhipiw 

wherein ^ corresponds to the convention term, K * corresponds to 

1 

Re 



the pressure gradient term, and corresponds to the diffusion term . 

15. (Currently Amended) A computer readable medium comprising a program 
stored thereon for numerical analysis of a flow field of incompressible viscous fluid, directly 
using V-CAD data, wherein the program causes causi n g a computer to perform the steps of : 
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a dividing step (A) e£dividing external data into a plurality of cells-(4-3} having 
boundaries orthogonal to each other, the external data including boundary data of an object 
that which contacts incompressible viscous fluid; 

a cell classifying step (B) e£classifying the divided cells into an internal cell (13a) 
positioned inside or outside the object and a boundary cell (13b) including the boundary data; 

a cut point determining step (C) ef-determining cut points in ridges of the boundary 
cell on the basis of the boundary data; 

a boundary face determining step (D) ef-determining a polygon connecting the cut 
points to be cell internal data for the boundary face; and 

a analyzing step (E) of applying a cut cell finite volume method combined with a 
VOF method to a boundary of a flow field to analyze the flow field , wherein step (E) 
comprises the steps of 

i. applying a two-dimensional QUICK interpolation scheme to a convection 
term for space integral; 

ii. applying a central difference having precision of a degree of a second 
order to a diffusion term; 

iii. combining the convection term and the diffusion term, and applying an 
Adams-Bashforth method having precision of a degree of a second order to the 
combined convection term and diffusion term for time marching; and 

iv. applying a Euler implicit method having precision of a degree of a first 
order to a pressure gradient term for time marching, 
wherein for a two-dimensional boundary cell, a governing equation in the finite 
volume method is expressed by a governing equation (7), 
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jj^dV = - jjdiv(u ® u)dV - jjdiv(pI)dV + — jjdiv(grad(u))dV 

/ 3/ ,/ Re y 



(7) 



"ljdiv{u0uW 

wherein corresponds to the convention term, 

^ corresponds to the pressure gradient term, and 



corresponds to the diffusion term . 
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